home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zlapll.z / zlapll
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. ZZZZLLLLAAAAPPPPLLLLLLLL((((3333FFFF))))                                                          ZZZZLLLLAAAAPPPPLLLLLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZLAPLL - two column vectors X and Y, let   A = ( X Y )
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE ZLAPLL( N, X, INCX, Y, INCY, SSMIN )
  13.  
  14.          INTEGER        INCX, INCY, N
  15.  
  16.          DOUBLE         PRECISION SSMIN
  17.  
  18.          COMPLEX*16     X( * ), Y( * )
  19.  
  20. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  21.      Given two column vectors X and Y, let
  22.  
  23.      The subroutine first computes the QR factorization of A = Q*R, and then
  24.      computes the SVD of the 2-by-2 upper triangular matrix R.  The smaller
  25.      singular value of R is returned in SSMIN, which is used as the
  26.      measurement of the linear dependency of the vectors X and Y.
  27.  
  28.  
  29. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  30.      N       (input) INTEGER
  31.              The length of the vectors X and Y.
  32.  
  33.      X       (input/output) COMPLEX*16 array, dimension (1+(N-1)*INCX)
  34.              On entry, X contains the N-vector X.  On exit, X is overwritten.
  35.  
  36.      INCX    (input) INTEGER
  37.              The increment between successive elements of X. INCX > 0.
  38.  
  39.      Y       (input/output) COMPLEX*16 array, dimension (1+(N-1)*INCY)
  40.              On entry, Y contains the N-vector Y.  On exit, Y is overwritten.
  41.  
  42.      INCY    (input) INTEGER
  43.              The increment between successive elements of Y. INCY > 0.
  44.  
  45.      SSMIN   (output) DOUBLE PRECISION
  46.              The smallest singular value of the N-by-2 matrix A = ( X Y ).
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.